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The funding provided in this grant supported the installa- 
tion of a Fabry-Perot Interferometer at the site of Arequipa, 
Peru. Measurements were made throughout the year to assess 
thermospheric winds and temperatures for an equatorial location. 
Appendix A reports on the analysis of the temperature data, with 
a comparison of the MSIS model for quiet and active magnetic 
periods, and for low and moderate levels of solar activity. 
Appendix B presents the results of thermospheric wind observa- 
tions throughout the same period. The chosen site of Arequipa, 
Peru is well suited to these observations because of the excel- 
lent seeing conditions, and the availability of high quality 
engineering support provided by the NASA lidar facility staff 
at this site. 

Although it was hoped that thermosperic winds and tempera- 
ture measurements could be made during the period of Project 
Condor in early March, 1983, this was not possible. The instru- 
ment was perational at Arequipa at the time requested, but cloud 
cover conditions prohibited any useful measurements. Some clear 
skies became available by mid-March, however, the full moon 
conditions during this time rendered the observations useless 
because of the high background associated with the lunar illumi- 
nation. The best measurements came during the dark moon periods 
of the following months beginning with April, 1983. Measurements 
were terminated in mid-August due to equipment problems caused 
by damage of the power invertor required to convert 50 hertz to 
60 hertz frequency. These problems have been remedied and the 
observatory is currently in operation at this time. 



Appendix A has been submitted to the Geophysical Research 
Journal for publication. The figures in the Appendix B were 
presented at the Thermospheric Dynamics workshop meeting at 
Calverton, Maryland, on the 4th of October, 1984. This material 
will form the basis for a publication in GRL. 
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MEASURED RESPONSE OF THE EQUATORIAL THERMOSPHERIC TEMPERATURE 
TO GEOMAGNETIC ACTIVITY AND SOLAR FLUX CHANGES 


Manfred A. Biondi 
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and 

John W. Meriwether, Or. 
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Abstract - Fabry-Perot determinations of thermospheric temperatures from 
630.0 nm nightglow line width measurements have been carried out for the 
period April-August 1983 from Arequipa, Peru (16.5°S, 71.5°W geographic; 

5°S magnetic). The nightly variation of the thermospheric temperature 
T n measured on 53 nights is compared with MSIS model predictions and 
found to agree occasionally with the model but, on average, to exceed 
model predictions by 180 K. The largest differences, 400- 500 K, 
occur during strongly increasing geomagnetic activity such as sudden 
commencements. Significant differences occur both during high geomagnetic/ 
low solar activity and during low geomagnetic/high solar activity. A 
suggested origin of the rapid increase in T n which follows sudden increases 
in geomagnetic activity with £ 1 day delay is the energy carried to 
equatorial regions from high-latitude (auroral oval) heat sources at 
both poles by gravity waves and equatorward neutral winds. 
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Introduction 

The temperature of the neutral thermosphere has been monitored 
by various ground-based and satellite techniques which include direct 
determinations such as measurements of the doppler width of the 630.0 nm 
airglow line and indirect methods such as radar backscatter plasma 
measurements and satellite orbit decay observations (see, for example, 
Jacchia, 1977; Hedin et al , 1977; Thuillier et al , 1977; Hernandez, 1982; 
Sipler et al, 1983 and the references cited therein). The temperature is 
expected to change in response to changing energy inputs in the heat 
balance equation governing the neutral thermosphere, in particular, to 
changes in the fluxes of solar radiation and of precipitating particles 
which are stopped in the upper atmosphere. A number of semi -empirical 
models (-e.g., Jacchia, 1977; Hedin et al , 1977; Thuillier et al , 1977) 
which predict thermospheric temperatures on the basis of indices related 
to solar euv and high-latitude precipitating particle fluxes, e.g. F] q . 7 » 
K p °T A p» etc., appear to provide reasonable agreement with measured 
nightly mean temperatures at mid-latitude (Hernandez, 1982); however, 
the measured diurnal variations and responses during geomagnetic storm 
conditions often are rather different from the predicted behavior. 

Recently, an automated airglow observatory has been set up at 
Arequipa, Peru (16.5°S, 71.5°W geographic, 5°S geomagnetic) in a joint 
effort by the Universities of Pittsburgh and Michigan, with the coopera- 
tion of the personnel of the NASA Laser Satellite Tracking Station at 
Arequipa. The principal observatory instrument is a field-widened 100 mm 
aperture Fabry-Perot interferometer suitable for nightglow 630.0 nm 
doppler width and doppler shift determinations. The unusually clear 
weather at Arequipa permits nightglow measurements on many successive 
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nights (the limitation in data-taking to date has largely been due to 
apparatus problems engendered by power line transients and the pervasive, 
fine dust carried by daytime winds). Thus, it has been possible to obtain 
sequences of ^ 12 hour measurements of nightime temperatures to investigate 
the short-term response of the thermosphere to changing solar euv radiation 
and geomagnetic activity. 

The present paper details measurements of the neutral temperature 
in the equatorial thermosphere over Arequipa for the period 1 April - 
9 August 1983, during which some 90 nights of useful data were obtained, 

53 of which have been reduced, and compares them with changes in solar 
flux and geomagnetic activity indices on both short- (few hour) and 
longer-term (weeks) time scales. 

Apparatus 

The optical apparatus and its electronic control hardware are 
very similar to those described in detail by Si pier et al (1983), while 
the computer control and data acquisition system and its software are 
of the type described by Meriwether et al (1983); therefore only a brief 
description will be given here. The system provides for automatic 
observations of the 630.0 nm nightglow line by the Fabry-Perot interferometer 
and a companion tilting filter photometer. The common 1 ines-of-sight of 
the two instruments are directed to the desired point in the sky by a 
two-axis, stepping-motor-driven pointing head on the roof of the observatory 
building. 

The tilting filter photometer monitors the wavelength range 
628-631 nm in order to detect interfering radiation (e.g., from OH) and 
cloud effects on the 630.0 nm observations. 

The 100 mm aperture Fabry-Perot (F-P) is index-of-refraction tuned 
("pressure-scanned") by changing the density of argon gas between the F-P plates 
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with a stepping-motor-driven volume changer consisting of a piston and 
cylinder sealed by a rolling diaphragm. The F-P is field-widened by 
a multi -aperture exit plate which provides a gain of ^ 3.8 in luminosity 
over a central aperture exit plate instrument. The 630.0 nm line profiles 
are obtained by coherently summing a number of wavelength scans across 
the line until either a predetermined signal level has been achieved or 
a maximum time interval has elapsed. After this, the computer directs 
the lines-of-sight of the instruments to the next observing position. 

Radiation from a stabilized 632.8 nm helium-neon laser is 
periodically (at intervals ranging from 20 to 120 min) introduced to 
the F-P to provide a determination of the instrumental broadening function, 
which is accomplished by computer- fitting the 632.8 nm line profile 
with a 10-term Fourier series. This truncated series is convolved with 
a gaussian line shape (the assumed form of the nightglow emission) and 
then least-squares fitted to the measured 630.0 nm profiles to yield 
line center frequency, line width, integrated intensity and background 
(Wickwar et al , 1984). The estimated error in the neutral temperatures 
determined from the line widths varies with the 630.0 nm nightglow 
intensity, ranging from - 30 K to - 70 K in most cases. 

Resul ts 

A summary of the neutral temperatures determined from the 630.0 nm 
doppler widths is presented in Fig. 1, together with published values 
of the solar flux (Sa) and geomagnetic planetary three-hour-range indices 
(Kp) from the World Data Center A for Solar-Terrestrial Physics. The 
nightly (- 12 hour duration) determinations of T n are indicated by the 
broad lines, most of whose width results from differences in temperature 
noted in different observing directions, e.g., N, E, S, W. The trend 
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of the lines indicates a rough average of the change of temperature 
during the course of each night, some nights yielding a modest fall, 
as normally expected, others indicating a sharper fall or a fall 
followed by a rise. The short, narrow lines on the T n plots are the 
MSIS model predictions (A. E. Hedin, private communication, 1984), 
using measured values of the solar and geomagnetic indices for the 
period shown. 

In order to compare short-term response of T n to changes in 
geomagnetic activity, detailed measurements on selected nights are 
shown in Fig. 2, with the symbols N, E, S, W indicating the observing 
azimuths (elevation = 30° in all cases). The solid lines are the 
MSIS model predictions. The sequence of six nights shown in Fig. 2 
involves a period of smoothly falling solar flux but a change from- 
geomagnetically quiet to disturbed back to quiet conditions. In passing, 
we note in the data of 4, 5, 6 August evidence of the post-midnight 
rise in T n which has been observed at low and equatorial latitudes by 
a number of investigators, for example, by Spencer et al (1979). 

Discussion and Conclusions 

It is clear from Figs. 1 and 2 that, while there are some nights 
when the T n values deduced from the 630.0 nm line width measurements 
agree reasonably with the MSIS model predictions, e.g., during August, 
many of the measurements yield temperatures significantly higher (100 K - 
400 K) than the predicted values. Non-thermal broadening of the 630.0 nm 
nightglow line by sources such as remnants of the dissociation kinetic 
energy from the electron-ion recombination which produces the 0(^D) atoms 
or kinetic energy of precipitating energetic oxygen atoms (Torr and Torr, 
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1984) can not account for the elevated temperatures measured. In the 
former case, complete slowing by collisions before radiation will occur 
at expected emission altitudes {< 400 km), while, in the latter, the 
relatively few energetic 0 atoms would have to share their kinetic energy 
with and simultaneously excite the ambient 0 atoms to the state in 
order to have a discernible effect near the doppler core of the line. 

We have been unable to find any instrumental effect which can 
account for erroneously high temperatures of the required magnitude. 

For example, spurious broadening due to instrumental drifts during the 
long periods of scan summing required by weak airglow emission can not 
account for the effect, since high temperatures were measured during 
periods of both low and high 630.0 nm nightglow intensity. Also, errors 
in instrumental width determinations due to improper frequency locking 
of the stabilized He-Ne laser would lead to an underestimate, rather 
than an overestimate of the 630.0 nm line width and, therefore, of T n . 

Significant differences between measured T n values and MSIS 
predictions occur both during periods of high solar/low geomagnetic 
activity, e.g., June 3-7 and during periods of low solar/high geomagnetic 
activity, e.g., April 1-16. Thus, the underestimate of the measured 
exospheric temperatures by the MSIS model can not be traced to an 
incorrect dependence on one activity index (solar or geomagnetic) alone. 

Our measurements of equatorial thermospheric temperatures over 
Kwajalein Atoll during August and September 1977 (Sipler et al , 1983) also 
yielded higher values, on the average, (by <v 330 K) than the MSIS model 
predictions. In addition, Hernandez (1982) has analyzed the T p values he 
obtained from 630.0 nm nightglow line width measurements at midlatitude 
(Fritz Peak Observatory, Colorado) from 1972 to 1979 and finds that his 
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measured values exceed the MSIS model predictions by an average of 100-150 K. 
The measured values in Fig. 1 exceed the MSIS model predictions by ^ 180 K 
on the average, so that substantial differences between values of T n 
determined by ground-based FP's and those of the MSIS model (empirically 
derived by fitting satellite mass spectrometer and incoherent scatter 
data to solar and geomagnetic indices) seem to be a continuing problem 
at mid- and equatorial latitudes. These differences can not be traced 
to an inherent, large discrepancy between ground-based F-P and incoherent 
scatter determinations, since simultaneous F-P and backscatter deter- 
minations of T n reported in two papers (Cogger et al , 1970 and Hernandez 
et al , 1975) revealed only a small difference, £ 30 K, for the T n 's 
determined by the two techniques. 

In examining the short-term response (£ 1 day) of the thermospheric 
temperature to rather rapid changes in solar EUV flux or to energy inputs 
at high-latitudes from particle precipitation, etc., it is inappropriate 
to use semi-empirical models such as MSIS or those of Jacchia, which 
are intended to describe average or longer-term behavior. Indeed, our 
findings of rapid responses of T n to strongly increased geomagnetic 
activity, especially sudden commencements, are not reflected in the MSIS 
model predictions for those periods (see, for example, the data of 
2-3 August, 6-8 August, and possibly 12-13 June in Figs. 1 and 2). In 
these cases, when a sudden increase in Kp occurs, the measured T n values 
remain low for up to a day and then increase sharply. 

A partial explanation of these observations may be found in 
gravity wave propagation of energy from a high latitude heating source 
to low latitudes. For example, recent calculations by Mayr et al (1984) 
have indicated substantial elevation of exospheric temperatures at 
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equatorial latitudes, AT n /T no ~ 0.2, several hours after the onset of 
Joule heating in the northern auroral oval (in the region of the heat 
source, C A W^no^ ma x ^ ”*)• Such efficient and rapid propagation of 
energy may well be the source of the occasional, rapid response of the 
low-latitude thermosphere to surges in high latitude heating. Also, 
equatorward winds driven by Joule heating in both the North and South 
auroral ovals and averaging 100-200 m/s would cause compressional 
heating of the equatorial thermosphere about a day after the onset. 

Fejer et al (1983) have attributed observed changes in the vertical 
drifts of the equatorial F-region over Jicamarca, Peru, which occur 
about one day after the onset of some geomagnetic storms, to similar 
changes in the thermospheric circulation pattern. These effects should 
be calculable by refinements of thermospheric general circulation 
models such as those of Dickinson et al (1981). 

In summary, the present measurements at equatorial latitudes 
yield average nighttime thermospheric temperatures which are signifi- 
cantly higher than MSIS model predictions, a finding similar to that 
noted at midlatitude by Hernandez (1982). Further, the rapid (£ 1 day) 
and large response of T n to sudden changes in geomagnetic indices 
suggests efficient and rapid transport of energy injected at high latitudes 
to the equatorial thermosphere, probably by gravity waves or by extensive, 
storm-generated equatorward winds. 
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Figure Captions 

Fig. 1 Nighttime thermospheric temperatures T n , geomagnetic 
planetary three-hour-range indices Kp and solar flux values Sa as 
a function of world date for March-August 1983. On the T n graphs: 
broad curves - measured values in various directions; short, narrow 
lines - MSIS model predictions. The arrows on the K p graphs 
indicate sudden commencements. 

Fig. 2 Nighttime thermospheric temperatures during August, 1983 

(OO^ 1 LT = 05^ UT). The three-hour K p indices are given at the top 

of each graph. Measured temperatures in various directions: 

+ + 

N, E, S, W symbols (the error bars, ranging from - 30 K to - 70 K, 
are omitted for clarity). MSIS model predictions: solid curves. 
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Arequipa FPI observatory: facts 


Field -"widened 10 cm FPI 
Automatic operations 

Located at Arequipa, Peru (16. 5S, 71.5W) 

(at the NASA lidar satellite tracking station) 

4.5 degrees south of magnetic equator 

FPI observes winds overhead of Jicamarca radar 

Observing directions typically 
Z.N.S.W or N.S.E.W 

Data selected from observations between 
new moon, April, 1983, and new moon, 

August, 1983; 55 nights reduced. 

Some results from 1984 but not all reduced. 

Error bars typically 15-20 m/s for winds, 

50 to 75 degrees for temperature 



Equatorial kinetic temperature: Results 

* Enhancement of FPI temperatures 
above quiet levels a few hours after 
the start of magnetic activity. 

Magnitude of this enhancement about 
300 degrees. This is followed by 
relaxation to pre-storm levels. 

* Apparent average offset of FPI 
temperatures from MS1S by 100-200 
degrees for relatively quiet times. 

* Definite suggestion of a midnight 
thermal enhancement for April and August 
data. Magnitude about 100 degrees. 

Seen in both 1983 and 1984 observatons. 



Hernandez: Mid-Latitude Thermospheric Neutjial Temperatures 



MS1S empirical model ( Hedin el ul.. I977|. 



exremuEWTii. temperature ik) 

Fig. II. Fit of the experimental data to four parameters.-.solar 
radio flux, geomagnetic activity, annual variation of the solar 
declination and a semi-annual variation. The correlation coefficient 
is 0.91. 



Fig. 4. Comparison of the experimental temperature with the 
empirical model of Thuillier et al. ( 1977a, fc|. 



EXPERIMENTAL TEMPERATURE IK) 

Fig. 5. Comparison of the experimental temperatures with the 
empirical model of Juccliia { 1977). 






INTENSITY TEMPERATURE (°K) ZONAL WINDS (m/s) 


300 

150 

0 

-150 

-300 


1800 

1500 

1200 

900 

600 


2.8 

2.2 

1.6 

1.0 


0.4 


3/4 JUNE 1983 4/5 JUNE 1983 6/7 JUNE 1983 




® m ®* 0 ® ® ® ® ®® 


. ®©«® ® ® 

© ® 

© ® 


18 21 24 3 618 21 24 3 618 21 24 3 6 

LOCAL TIME 





MERIDIONAL 

INTENSITY TEMPERATURE (°K) WINDS (m/s) 


16/17 APRIL 1984 19/20 APRIL 1984 



LOCAL TIME 








TEMPERATURES (°K) 


1800 
1200 
6 00 




LOCAL TIME 






INTENSITY 



15/16 APRIL 1983 5/6 JUNE 1983 23^24 JUNE 1983 30/31 JULY 1983 


MERIDIONAL 

INTENSITY TEMPERATURE (°K) WINDS (m/s) 


16/17 APRIL 1984 19/20 APRIL 1984 



LOCAL TIME 







ZONAL WINDS (m/s) 









MERIDIONAL WINDS (m/s) 


300 

150 

0 

-150 


t — i — | — i — i — i — i — i — i — r 
31 JULY/I AUG. 1983 


4~i~“ 

k i 

J I I I I 1 ! I I 1 1 


•300 

300 1 — i — i — t - ! — i — i — i — i — r 

1 1 3/4 AUG. 1983 


150 

0 

-150 

-300 

300 

150 

0 

-150 

-300 


1 T 


__ J 

-<j> <p 


I . 1 


J I I L 


1 5/6 aIjg! 198^ 


i — r 


- ----- ' ^---4— ■ 




I I I I I I I 1 I L 


"i i | i i — i — i — i — 7 — r 
2/3 AUG. 1983 


■f f 


J 1 l_J I 1 L 


J I L 


t — I — I — r 


1 I I I 

4/5 AUG. 1983 


r— i — Zi — i — r 






j i i i 


J I I L 


i — i — | — i — i — r~i — i — i — i — r 
6/7 AUG. 1983 





j i l i i t i i l „ . i 


18. 21 24 3 6 18 21 24 3 

LOCAL TIME 












ZONAL WINDS (m/s) 


300 


150 

0 

-150 

-300 



18 21 24 3 6 18 21 24 3 6 


LOCAL TIME 











ZONAL WINDS (m/s) 










MERIDIONAL 

INTENSITY TEMPERATURE (°K) WINDS (m/s) 


16/17 APRIL 1984 



LOCAL 






MERIDIONAL WINDS (m/s) 


300 

150 


0 

-150 


-300 



18 21 24 3 6 18 21 24 3.6 

LOCAL TIME 







ZONAL WINDS (m/s) 



LOCAL TIME 







INTENSITY TEMPERATURE (°K) ZONAL WINDS (m/s) 


300 

150 

0 

-150 

-300 


3/4 JUNE 1983 4/5 JUNE 1983 6/7 JUNE 1983 


* t 


1800 

1500 

1200 

900 

600 


2.8 

2.2 


0.4. 


® o® ® ® ® 



21 24 3 618 21 24 3 618 21 24 3 6 

LOCAL TIME 








MERIDIONAL WINDS (m/s) 




LOCAL TIME 











MERIDIONAL WINDS (m/s) 



LOCAL TIME 






Intensity Temperature (k) Ulnd (m/s) 


Horizontal wind, temperature, and intensity observations 
DAY 271, 1984 
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CONCLUSIONS 


6300A surface brightnesses: Results 

* Definite enhancement of 6300A 
surface brightness in the south as 
compared with other directions. 
Probably connected to the tropical 
airglow arcs. Meridional winds 
small (<25 m/s) throughout night. 

* Indication of northward migration 
of the observed 6300A enhancement 
in the evening hours as 
observations approach local winter 
solstice. This is probably related to 
the observed poleward (to the south) 
meridional wind (of magnitude 50 m/s) 
in this period. 



Equatorial thermospheric winds: Results 

* Zonal component of winds always 
eastward, but speed approaches zero 
sooner near equinox than at summer, 
solstice. Typical magnitude at peak 

is of order 100 to 150 m/s. 

* Suggestion of zonal wind increase 
after twilight and recovery of 6300A 
emission for April data. Origin not 
clear but may be related to mid-night 
thermal enhancement. 

* Meridional wind virtually zero for 
equinox in 1983; shows evening 
flow towards winter hemisphere in 
early evening for solstice data. 
Suggestion of post-midnight surge in 
April 1984 data. 

* No major effects associated with 
mag. storm activity. Suggestion of 
decrease in zonal component below 
nominal levels. 



